Wrinkling hierarchy in constrained thin sheets from suspended graphene to curtains 
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We show that thin sheets under boundary confinement spontaneously generate a universal self- 
similar hierarchy of wrinkles. From simple geometry arguments and energy scalings, we develop a 
formalism based on wrinklons, the transition zone in the merging of two wrinkles, as building-blocks 
of the global pattern. Contrary to the case of crumple paper where elastic energy is focused, this 
transition is described as smooth in agreement with a recent numerical work [1 . This formalism 
is validated from hundreds of nm for graphene sheets to meters for ordinary curtains, which shows 
the universality of our description. We finally describe the effect of an external tension to the 
distribution of the wrinkles. 



The drive towards miniaturization in technology is 
demanding for increasingly thinner components, raising 
new mechanical challenges [2 . Thin films are however 
unstable to boundary or substrate-induced compressive 
loads: moderate compression results in regular wrin- 
kling [3i7j while further confinement can lead to crum- 
pling [8l|9]. Regions of stress focusing can be a hindrance, 
acting as nucleation points for mechanical failure. Con- 
versely, these deformations can be exploited construc- 
tively for tunable thin structures. For example, singular 
points of deformation dramatically affect the electronic 
properties of graphene [10 . 

Here, we show that thin sheets under boundary con- 
finement spontaneously generate a universal self-similar 
hierarchy of wrinkles; from strained suspended graphene 
to ordinary hanging curtains. We develop a formalism 
based on wrinklons, a localized transition zone in the 
merging of two wrinkles, as building-blocks to describe 
these wrinkled patterns. 

To illustrate this hierarchical pattern, we show in 
Fig. [l^ wrinkled graphene sheet along with an ordinary 
hanged curtain. These patterns are also similar to the 
self-similar circular patterns first reported by Argon et 
al. for the blistering of thin films adhering on a thick 
substrate [IT. The diversity and complexity of those 
systems, characterized by various chemical and physical 
conditions, could suggest, a priori^ that the underlying 
mechanisms governing the formation of these patterns 
are unrelated. However, these systems can be depicted, 
independently from the details of the experiments, as a 
thin sheet constrained at one edge while the others are 
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free to adapt their morphology. These constraints can 
take the form of an imposed wavelength at one edge or 
just the requirement that it should remain flat. 

As illustrated in Fig.[l]and[2j sheets made from various 
materials constrained at one edge by an imposed sinu- 
soidal profile spontaneously develop a hierarchical pat- 
tern of folds or wrinkles. At first sight, as quoted by 
numerous authors [H fTTHTT] , these patterns consist of a 
hierarchy of successive generations of folds whose typical 
size gradually increases along x (Fig. ^p)- We propose 
to rationalize these various hierarchical patterns by con- 
sidering the evolution of the average wavelength. A, with 
the distance to the constrained edge, x. This evolution 
is adequately described by a simple power law, A ~ x^, 
see Fig. [l]^, which confirms the self-similarity of these 
patterns as hypothesized in previous theoretical studies 
j9| [T2HT4] . Interestingly, curtains made of various materi- 
als with contrasted properties exhibit similar exponents. 
We observe values close to 2/3 for "light" sheets and to 
1/2 for "heavy" sheets. Therefore the exponent m is a 
robust feature of these folding patterns. 

In this work, we describe in terms of simple scaling 
laws the theoretical arguments developed in the mathe- 
matical studies of the von Karman equation [9l [T3l [T4l [T6] 
and infer the properties of the hierarchical patterns. We 
also compare these results with extensive experimental 
data. To the best of our knowledge, the experimental 
characterization of these patterns had not been carried 
out before. 

Assuming inextensibility of the sheet along the y direc- 
tion, the imposed undulation along this direction exactly 
compensates for an effective lateral compression of the 
membrane by a factor (1 — A) defined as, (1 — A) = 

W/Wq = W/ y^l^ \dz/dy)^dy, where Wo and W 
are the curvilinear and projected width of the curtain, re- 
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FIG. 1: a, SEM image of a graphene bilayer thin sheet 
suspended across pre-defined trenches on Si/Si02 substrates 
(scale bar: l/xm). b, Hierarchical pattern of folds obtained 
for a long suspended curtain made of a thin sheet of rubber 
(scale bar: 25 cm), c, Evolution of the average wavelength, 
A, with the distance from the constrained edge, x, for vari- 
ous curtains as indicated in the legend. Power law fits are 
added (the power exponents, m, are close to 2/3 for the short 
fabric and the paper curtains and 1/2 for the long fabric and 
the rubber curtains). Inset: Evolution of A with x for the 
graphene sheet (m is equal to 0.45 ± 0.02). The experimen- 
tal parameters are detailed in the supplementary information. 



spectively, and z{x^y) is the out-of-plane deformation of 
the sheet. At any position along the x axis, the function 
z{x^y) is typically sinusoidal along with an amplitude 
A{x) and a wavelength X{x). The inextensibility hypoth- 
esis along the y axis imposes A ~ (A/A)^ at the lowest 
order, where the lateral compression is assumed to be 
constant throughout the length of the curtain. The un- 
dulations of the sheet along y are characterized by a cur- 
vature hi d'^z/dy'^ whose typical value, varying along 
X, is of order hz{x) ~ A/X^. The corresponding energy 
per unit area, iib, for bending the membrane is thus of 
order ui^ ^ Eh^ ti^ ^ Eh^A/X^, where E is the Young 
modulus and h the thickness of the sheet. Since u\y is 
proportional to 1/A^, the membrane adopts the largest 
possible wavelengths, in order to minimize energy. This 
tendency to increase the wavelength, combined with the 
constraint imposed at the boundaries, is the source of the 
observed hierarchical wrinkling pattern. 

Inspired by previous models based on successive 
period-doubling transitions [H [H [M], we consider that 
the allometric laws mentioned above can be derived by 
considering that the global pattern results from the self- 
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FIG. 2: a. Schematic representation of the wrinklon exper- 
iments, b. Morphology of the transition A to 2A for a con- 
strained plastic sheet for A = 6 mm and A = 8 mm. c. 
Evolution of the normalized length of a wrinklon, L/X, with 
the normalized amplitude, A/h (fixed wavelength, A = 8 mm) 
for different thicknesses as indicated. 



assembly of building-blocks which we denote wrinklons. 
A single wrinklon corresponds to the localized transition 
zone needed for merging two wrinkles of wavelength A 
into a larger one of width 2A. This transition requires 
a distortion of the membrane which relaxes over a dis- 
tance L. In other words, each wrinklon is characterized 
by a size, L, which depends on the material properties 
and on the wavelength A. To investigate the proper- 
ties and behavior of wrinklons, we have performed model 
experiments using thin plastic sheets. The sheets were 
constrained with sinusoidal clamps: two opposite edges 
are constrained by a wavelength A (amplitude A) and 2A 
(amplitude 2A), respectively, see Fig. [2|L,b. The normal- 
ized size of the wrinklons, L/A, is plotted in Fig. ^ as 
a function of the normalized amplitude, A/h; the data 
collapse on a single curve defined by L/A - y^A/h. This 
relation implies that L ex A^/^ since A ~ X^/A. 

In a further step, the wrinklons can be assembled to 
mimic the behavior of a complete hierarchy. Indeed, if L 
is the distance over which the wavelength increases from 
A to 2 A, its variation, dA/dx, is thus of order X/L. Hence, 
the evolution of A as a function of the distance from the 
constrained edge, x, is given by 



dA 

dx 



X 
L' 



(1) 



Considering the scaling L cx A^/^ deduced from the sin- 
gle wrinklon experiments, equation ([T]) indicates that the 
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wavelength along the sheet should evolve like A oc x^l'^ . 
The excellent agreement between this power law and the 
experimental data measured for light sheets (Fig.jl]^) pro- 
vides a strong support to the concept of wrinklons as 
building-blocks. Equation ([T]) can now be regarded as a 
tool that connects the properties of single wrinklons to 
the features of the full wrinkling-cascade pattern. 

We now focus on the description of an elementary 
building block. For confined thin sheets, stretching de- 
formations are costly as compared to pure bending. The 
sheet tends to adopt an isometric (developable) shape [8] . 
However, the only developable solutions compatible with 
boundary conditions generally include flat domains sur- 
rounded by edge or point-like singularities. These sin- 
gularities, which focus the elastic energy into narrow re- 
gions, have been classified as developable cones [18|fT9], 
ridges [3 [20], or curved ridges [21 . In our case, the 
scenario is however significantly different: in contrast to 
crumpling, stretching is smoothly distributed in the tran- 
sition zone as pointed out recently in numerical simula- 
tions of deformed membranes [1 . The necessary stretch- 
ing required for connecting the periodic patterns can 
be illustrated by a simple origami model made with a 
sheet of paper (see supp. info.). The stretching energy 
can be estimated through the elongation strain of the 
sheet along x within a transition domain. The typical 
value of the strain along x is of order a^, where a ~ 
AjL X/S}^'^ /L is the average slope of the membrane. 
The stretching energy thus reads ^ Eh{a^)'^L\ ^ 

As observed in Figs, [l] and [2j wrinklons should also 
include a tip singularity (a small region where Gaus- 
sian curvature is large). This singularity can be de- 
scribed as a semi-circular fold of radius p (Fig. ^p)- 
The energy of these singularities has been derived by 
Pogorelov [21] in a study of deformed shells. In our 
context, the energy of such curved folds reads ^ 
Eh^/^a^/^p^/^ - Eh^/^ A^/^X^^L-^, where the radius 
at the tip of the wrinklon is taken as p ~ A^/L as sug- 
gested by the roughly parabolic shape of the crest of the 
defect (Fig. ^p)- Nevertheless, the ratio of the curved 
fold energy to the stretching energy of the wrinklon, 
Uci/Us ~ (/i/A)^/^, is very small in our experiments : 
the effect of this concentrated region can therefore be 
neglected in the following. 

The total energy of a wrinklon, of characteristic area 
LA, is thus given by Utot = t/s + - EhX^/S^L'^ + 
Eh^ALX~^. The size of a single wrinklon is finally ob- 
tained by minimizing /7tot with respect to L, yielding 



L(A) ~ 



(2) 



This scaling emerges from a balance between bending 
and stretching energies and was previously reported for 
other situations, such as the decay length of an imposed 
curvature in a sheet [20^ or the extension of a pinch in 
a pipe [22]. The scaling for the wavelength describing 
the whole hierarchical pattern is obtained by integration 
of equation ([T]) with L{X) given by equation ([2| and is 
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FIG. 3: Master curves gathering all data, a, Normalized 
wavelength, A = XA^^^/h, as a function of the normalized 
distance, x = x/h, from the constrained edge for short light 
sheets (fabric curtain, paper curtain and constrained plastic 
sheets). Dashed line: A = 2.89 b, Normalized wave- 

length A = A//i as a function of the normalized distance from 
the constrained edge x = {x/h){Eh/Ty^'^ for sheets under 
tension: fabric curtains, rubber curtains, suspended bilayer 
graphene sheet and polystyrene thin films deposited on water 
from Ref. [H]. Dashed line: A = 2.85 x^'^^. 



found to be 



A(x)AV^ 
h 



2/3 



(3) 



The scaling law, A ex x^/^, is in very good agreement with 
the observed power laws for light curtains, e.g. made of 
fabric or paper sheets (Fig. [l]^). In addition to yielding 
the proper exponent, this relation enables the comparison 
of the data obtained from seemingly disparate systems, 
over a wide range of lengthscales and independently of 
material properties. Figure provides a remarkable col- 
lapse of the evolutions of the wavelengths measured with 
light curtains and various thin plastic sheets. 

Heavy curtains, made from fabric or rubber, and 
constrained graphene bilayers do not follow this be- 
haviour (instead, they obey A ex x^/^). In these ex- 
periments, an additional tensile force is acting on the 
sheet. This tension, T, is given by the longitudinal ten- 
sile strain induced by thermal manipulation in the case 
of graphene sheets [6 and by gravity for heavy curtains 
(T = pcgh{H — x) ^ pcghH, where pc, 9^ and H 
are the density of the curtain, the gravity constant, the 
thickness and the height of the curtain). These systems 
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can also be compared to the cascade of wrinkles observed 
for compressed thin polystyrene films on an air-water in- 
terface [15] since the surface tension of water at the free 
edges pulls the thin sheet. 

The tension exerted along x per unit width, imposes 
an additional stretching energy given hy Ut ^ T o? LX ~ 
TAA^L"^, and becomes dominant when Ut > /7s, that 
is when T > Eh'^A/A. The total energy of the distorted 
membrane thus becomes [/tot = Ut -\-U\). The length of 
a wrinklon found from the minimization of [/tot is 



(4) 



Similar relations, reflecting a balance between tension 
and bending energies, were previously proposed for sin- 
gle wavelength patterns in stretched sheets and heavy 
curtains [11 [17]. As expected, the tensile force increases 
the length of wrinklons for a given wavelength [16]. By 
integration of equation ([T]) with L{X) given by equa- 
tion Q, we obtain the corresponding spatial evolution 
of the wavelength along a heavy sheet 
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(D 
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(5) 



This scaling is in excellent agreement with the power 
laws observed for heavy curtains and graphene bilay- 
ers, (Fig. [l]^). The data of various macroscopic cur- 
tains, graphene bilayers and nanometric polystyrene in- 
deed collapse onto a single master curve without any fit- 
ting parameters (see Fig. Our formalism is thus 
validated from hundreds of nm for graphene sheets to 
meters for rubber and fabric curtains, which shows the 
universality of our description. The transition between 
the stretching and tension regimes can be obtained by 
comparing the relations (|3| and ([5|. The critical dis- 
tance from the edge at which this transition occurs is 
given by x^/h ~ {Eh/T)^^'^A. In gravity dominated 
systems, the tension T :^ pghH, gives the typical cur- 
tain length He ~ h{E/ pgh)^^^A'^^^ above which tension 



dominates. Curtains shorter than He (about Im for our 
fabric) were used to observe the regimes dominated by 
stretching ("light sheets"), whereas the top part of longer 
curtains were used for experiments concerning "heavy 
sheets" . 

In summary, we showed that the self-similar patterns 
observed in sheets constrained at one edge cannot be de- 
scribed with d-cone or ridges singularities. In contrast, 
they can be built by stitching together building-blocks, 
called wrinklons characterized by a diffuse stretching en- 
ergy. The self-similar structure is then related to the size 
of these wrinklons that depends on material properties 
and the local wavelength. Interestingly, we also show that 
these building-blocks can be readily manipulated through 
the size and energy cost of a single wrinklon by applying 
a tension. For large values of tension, we even expect a 
transition towards a purely cylindrical pattern along the 
sheet with a single wavelength. Finally, we can draw a 
parallel between this study and the previously reported 
fractal buckling of torn plastic sheets [27. The imposed 
metric indeed determines the three-dimensional shape of 
the distorted membrane, characterized by a superimposi- 
tion of various modes. In contrast, the patterns observed 
here for constrained thin sheets exhibit a continuous evo- 
lution of the wavelength. 
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Supplementary information for "Wrinkling hierar- 
chy in constrained thin sheets from suspended 
graphene to curtains" 

I. METHODS 

A. Curtains 

The curtain experiments were carried out using fabric: 
thickness 220 ^m, density pc = 820 kg/m^ and elastic 
modulus E ^ 1 MPa; natural latex: thicknesses 220 and 
400 /im, density pc = 980 kg/m^ and elastic modulus 
E ^ 1 MPa; paper: thickness 125 pm. The fabric cur- 
tain was 4 m long and 2.5 m wide, while the rubber 
and paper curtains were 2 m long and 1 m wide. Each 
curtain was hung with an imposed sinusoidal boundary 
condition through an array of screws of tunable length. 
The wavelength was fixed to A = 20 mm with an ampli- 
tude A spanning from 1.2 to 8 mm. The bottom edges 
of the curtains were left free except for the ballasted cur- 
tain experiments where 10 steel disks of 180 g each were 
hung from the bottom edge with a spacing of 100 mm 
in between. The average wavelength was obtained by di- 
viding the total width of the curtain by the number of 
wrinkles found at a given distance from the constrained 
edge. For the various curtains, the evolution of A with 
the distance from constrained edge was analyzed with 
power laws. These fits yield the following exponents: 
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wafers with pre-patterned trenches. Bilayer graphene 
sheets were identified by color contrast in an optical mi- 
croscope and/or Raman spectroscopy described in W. 
Bao et al Nature Nanotech. 4, 562 (2009). The Si 
substrates were p-doped, and the thickness of silicon and 
Si02 are 0.5 mm and 300 nm, respectively. The trenches, 
fabricated at the UCSB Nanofabrication facility, were de- 
fined by photolithography followed by plasma etching in 
a reactive ion etcher (RIE) system. Strain in graphene 
sheets was estimated using the same procedure as the one 
used to produce Fig. le in W. Bao et al. (2009). 

SEM image of graphene bilayer thin sheets suspended 
across pre-defined trenches on Si/Si02 substrates and an- 
nealed at 523K were recorded with a tilt angle of 75° to 
improve the contrast in amplitude. 



C. Isolated wrinklons. 

These experiments were carried out using A4 size 
sheets made of biaxially oriented polypropylene (BOPP 
Innovia) of thicknesses 30, 50, 90, 100 and 250 pm 
and elastic modulus ^ ^ 2.6 GPa. To impose a sinu- 
soidal condition at the edge of the sheet we used Plex- 
iglas clamps of thickness 5 mm with sinusoidal profiles 
of different amplitudes and wavelengths to clamp the 
sheets: fixed amplitude A = ?> mm with wavelengths 
A = 13, 18, 26 and 58 mm; fixed wavelength A = 24 mm 
with amplitudes A = 3, 4, 5, 6 and 7 mm. A second 
clamp of wavelength 2A and amplitude 2A was used to 
impose a second profile condition for the sheet at a cer- 
tain distance L2 away from the clamped edge. In this 
way we forced the sheet to switch from the initial wave- 
length A to the next wavelength 2 A, thus generating the 
formation of a controlled wrinklon. The distance L2 was 
selected by separating progressively the opposite clamps 
until the shape of the wrinkles became stationary and 
before a new generation of wrinklons could arise. 



B. Graphene bilayers 



Suspended graphene membranes were prepared by 
standard mechanical cleavage technique on Si/Si02 



D. Morphology of wrinklons. 
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FIG. 4: a, Origami describing the distorted domain at the A to 2 A transition. This toy model is illustrative of a sheet with a 
vanishing thickness and shows the necessity to stretch the sheet to connect both A and 2 A patterns. For sheets with vanishing 
thickness, the deformation energy is concentrated in sharp singularities (ridges or d-cones). The width/size of these singularities 
decreases as the thickness tends to zero. The limit would then correspond to sharp folds that is modeled by simple origami. 
The scale bar is 5 cm. b, Schematic representation of the transition domain for sheets with finite thickness. The slope, a, along 
X direction is indicated. 



